arXiv: 1502.0415 lv2 [physics.optics] 17 May 2015 


Electromagnetically induced transparency for guided light in an atomic array outside 

an optical nanofiber 

Fam Le Kien 

Wolfgang Pauli Institute, Oskar Morgensternplatz 1, 1090 Vienna, Austria 

A. Rauschenbeutel 

Vienna Center for Quantum Science and Technology, Institute of Atomic and Subatomic Physics, 

Vienna University of Technology, Stadionallee 2, 1020 Vienna, Austria 
(Dated: May 19, 2015) 

We study the propagation of guided light along an array of three-level atoms in the vicinity of 
an optical nanofiber under the condition of electromagnetically induced transparency. We examine 
two schemes of atomic levels and field polarizations where the guided probe field is quasilinearly 
polarized along the major or minor principal axis, which is parallel or perpendicular, respectively, 
to the radial direction of the atomic position. Our numerical calculations indicate that 200 cesium 
atoms in a linear array with a length of 100 pm at a distance of 200 nm from the surface of a 
nanofiber with a radius of 250 nm can slow down the speed of guided probe light by a factor of 
about 3.5 x 10 6 (the corresponding group delay is about 1.17 ps). In the neighborhood of the Bragg 
resonance, a significant fraction of the guided probe light can be reflected back with a negative group 
delay. The reflectivity and the group delay of the reflected field do not depend on the propagation 
direction of the probe field. However, when the input guided light is quasilinearly polarized along the 
major principal axis, the transmittivity and the group delay of the transmitted field substantially 
depend on the propagation direction of the probe field. Under the Bragg resonance condition, an 
array of atoms prepared in an appropriate internal state can transmit guided light polarized along 
the major principal in one specific direction even in the limit of infinitely large atom numbers. The 
directionality of transmission of guided light through the array of atoms is a consequence of the 
existence of a longitudinal component of the guided light field as well as the ellipticity of both the 
field polarization and the atomic dipole vector. 

PACS numbers: 42.50.Gy, 42.50.Nn, 42.81.Dp, 42.81.Gs 


I. INTRODUCTION 

Optical properties of materials can be dramatically 
modified by quantum interference between the excita¬ 
tion pathways 0,i- Intensive attention has been de¬ 
voted to the manipulation and control of the propagation 
of light through coherently driven optical media, espe¬ 
cially in the connection with the possibility of enormous 
slowing down, storage, and retrieval of optical pulses [lj— 
H]. The interest to this topic is related to the applica¬ 
tions for optical delay lines, optical data storage, optical 
memories, quantum computing, and sensitive measure¬ 
ments. Through the technique of electromagnetically in¬ 
duced transparency (EIT) 0-0, ultralow group velocities 
of light have been obtained in hot [zj and cold 0 atomic 
gases. This technique allows one to render the material 
highly transparent and still retain the strong dispersion 
required for the generation of slow light. In addition, the 
transmitted pulse can experience strong nonlinear effects 
due to the constructive interference in the third-order 
susceptibility x^. This leads to new techniques in non¬ 
linear optics at the few-photon level, which may find im¬ 
portant applications to quantum information processing. 

The EIT technique has been extended to media embed¬ 
ded in a variety of waveguide systems, such as rectangular 
waveguides [l(|, hollow-core photonic-crystal fibers E3- 
Il3| . coupled resonator optical waveguides [14], quantum 


well waveguides M , waveguide-cavity systems pj| , and 
optical nanofibers EMI. EIT-based photon switches in 
waveguides [l3l[23l] have been examined. The generation 
and waveguiding of solitons in an EIT medium have been 
studied [24| . Waveguiding of ultraslow light in an atomic 
Bose-Einstein condensate [25| has been investigated. 

Nanofibers are optical fibers that are tapered to a di¬ 
ameter comparable to or smaller than the wavelength of 
light |26h28I ]. Slowing down of guided light in an opti¬ 
cal nanofiber embedded in an EIT medium has been in¬ 
vestigated jl7l - f22t | . The first observation of EIT at very 
low power levels of guided pump and probe light was 
reported in Ref. [20|. Very recently, coherent storage 
of g uided light has been experimentally demonstrated 
Hi HI- In Ref. El, the propagation of light is de¬ 
scribed in terms of the averages of the local refractive in¬ 
dex, the local absorption coefficient, and the local group 
delay in the fiber cross-section plane. The studies in 
Refs, ii are based on a more rigorous formalism 
that takes into account the inhomogeneous density dis¬ 
tribution of the atomic gas and the inhomogeneous mode- 
profile function of the guided field in the fiber trans¬ 
verse plane. The atomic medium considered in iT7l - l2l] is 
continuous. Meanwhile, recent experiments with atom- 
waveguide interfaces |29kUI used linear arrays of atoms 
prepared in a nanofiber-based optical dipole trap [[35l — 
1371 ]. It has recently been demonstrated experimentally 
that spin-orbit coupling of guided light can lead to direc- 
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tional spontaneous emission [dU [HU and optical diodes 
ld~H . When the array period is near to the Bragg reso¬ 
nance, the discreteness and periodicity of the array may 
lead to significant effects, such as nearly perfect atomic 
mirrors, photonic band gaps, long-range interaction, and 
self-ordering [39l - l5ll | . In the prior work on atoms trapped 
in a one-dimensional optical lattice under the EIT con¬ 
dition [HU HU , scalar light fields in free space were con¬ 
sidered. Scattering of a scalar light field from an array 
of three-level atoms with two degenerate lower levels in 
a waveguide has also been studied [52). 

In a nanofiber, the guided field penetrates an appre¬ 
ciable distance into the surrounding medium and appears 
as an evanescent wave carrying a significant fraction of 
the propagation power and having a complex polariza¬ 
tion pattern [HU HU- Since the nanofiber is thin, the 
guided modes of the nanofiber are the fundamental HEn 
modes [HU [HU- These modes are hybrid modes. The 
field in such a mode has longitudinal electric and mag¬ 
netic components. The local polarization of the mode 
varies in the fiber cross-section plane [HU and depends 
on the propagation direction [33l , 138| j. Therefore, the use 
of the scalar field formalism to treat the interaction of a 
nanofiber-guided field with an atom is not always appro¬ 
priate. When the local polarization of the guided field 
is elliptical and the dipole matrix-element vector of the 
atom is a complex vector, direction-dependent effects in 
the atom-field interaction may occur [HU [HU, IS, [HI, [HH] . 

In view of the recent results and insights, it is necessary 
to develop a systematic theory for the propagation of 
guided light under the EIT condition in an atomic array 
taking into account the vector nature of the guided field 
and the discreteness and periodicity of the array. 

In this paper, we study EIT in a one-dimensional peri¬ 
odic array of three-level atoms trapped along an optical 
nanofiber. We examine two schemes of atomic levels and 
field polarizations where the guided probe field is quasi- 
linearly polarized along the major principal axis x or the 
minor principal axis y , which lie in the fiber cross-section 
plane and are parallel or perpendicular, respectively, to 
the radial direction of the atomic position. We take into 
account the vector nature of the guided field and the dis¬ 
creteness and periodicity of the atomic array. We study 
the transmittivity and reflectivity of guided light and the 
time evolution of the transmitted and reflected fields. 

The paper is organized as follows. In Sec. [TT] we de¬ 
scribe two schemes of atomic levels and field polariza¬ 
tions for nanofiber-based EIT, and present the coupled¬ 
mode propagation equations. In Sec. lIIIl we study EIT in 
the homogeneous-medium approximation and the phase¬ 
matching approximation. In Sec. IIVI we investigate EIT 
in a discrete array of atoms with the help of the transfer 
matrix formalism. Our conclusions are given in Sec. [V] 


(a) x-polarized probe z q 00 
z and c-type transition -f 



(b) y-polarized probe z q 00 
z and 7r-type transition -f 
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FIG. 1. (Color online) Left column: An array of atoms out¬ 
side an optical nanofiber with a plane-wave control field £ c , 
a guided probe field £ p , and a guided reflected field £b . The 
atomic array is aligned parallel to the fiber axis z and lies 
in the zx plane of the Cartesian coordinate frame {x,y,z}. 
The probe field £ p is quasilinearly polarized along the ma¬ 
jor principal axis x in scheme (a) and the minor princi¬ 
pal axis y in scheme (b). The control field £ c propagates 
along the axis x and is linearly polarized along the axis y 
in scheme (a), and propagates along the axis y and is cir¬ 
cularly polarized in scheme (b). Right column: The dia¬ 
gram of atomic energy levels and electric dipole transitions 
for the analysis. Atomic cesium is used. The quantization 
axis zq is the axis y. In both schemes (a) and (b), the 
upper level is |e) = \ 6 P 3 / 2 ,F' = 4 , M' = ±4). In scheme 
(a), the lower levels are | g) = | 65 i/ 2 , F = 3, M = ±3} and 
| h) = 1651/2, F = 4, M = ±4). In scheme (b), the lower lev¬ 
els are \g) = \ 6 Si/ 2 ,F = 4, M = ±4) and | h) = \ 6 Si/ 2 ,F = 
3, M = ±3). 


II. NANOFIBER-BASED EIT SCHEMES AND 
COUPLED-MODE PROPAGATION EQUATIONS 

We consider a linear periodic array of A-type three- 
level atoms trapped outside an optical nanofiber (see Fig. 
[T]). The nanofiber has a cylindrical silica core, with radius 
a and refractive index m, surrounded by vacuum, with 
refractive index n 2 = 1. We use the Cartesian coordinate 
system { x , y , z} and the associated cylindrical coordinate 
system {r,ip,z}, with z being the fiber axis. We assume 
that the array of atoms is parallel to the fiber axis z 
and lies in the zx plane. The positions of the atoms in 
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the array are characterized by the Cartesian coordinates 
Xj = Xq, yj = 0, and Zj = (j — 1)A. Here, the index 
j = 1, 2,..., N labels the atoms, with N being the num¬ 
ber of atoms in the array, and the parameter A is the 
period of the array. The axes x and y , which lie in the 
fiber cross-section plane and are parallel and perpendic¬ 
ular, respectively, to the radial direction of the atomic 
position, are called the major and minor principal axes, 
respectively. Although our theory is general and applica¬ 
ble, in principle, to arbitrary multilevel atoms, we assume 
cesium atoms throughout this paper. For simplicity, we 
neglect the effect of the surface-induced potential on the 
atomic energy levels. This approximation is reasonable 
when the atoms are not close to the fiber surface. We 
also neglect the effect of the far-detuned trapping light 
fields. 


A. Quasilinearly polarized nanofiber-guided modes 

We represent the electric component of a nanofiber- 
guided light field as E = (£e~ lult + c.c.)/2 = ( £\ie~ lult + 
c.c.)/2, where u> is the angular frequency and £ = £u 
is the slowly varying envelope of the positive-frequency 
part, with £ and u being the field amplitude and the 
polarization vector, respectively. In general, the ampli¬ 
tude £ is a complex scalar and the polarization vector u 
is a complex unit vector. The guided light field can be 
decomposed into a superposition of quasilinearly polar¬ 
ized modes [53|. These guided modes can be labeled by 
the index (w,/, £), where / = +1 or —1 (or simply + or 
— ) stands for the positive (+z) or negative (— z) prop¬ 
agation direction, respectively, and £ = x or y stands 
for the major polarization axis. In the cylindrical co¬ 
ordinates, the transverse-plane profile functions of the 
positive-frequency parts of the electric components of the 
modes (w, /, £) are given by [HI, Hl-fHi] 

=v / 2 (re r cos p + ipe^ sin p + fze z cos <p), ^ 

e (^/y) =y /2 (re r sin ip — i(f>e v cos p + fze z sin ip). 

Here, the notations r = x cos p + y sin p, ip = —x sin p + 
y cos ip, and z stand for the unit basis vectors of the cylin¬ 
drical coordinate system, where x and y are the unit 
basis vectors of the Cartesian coordinate system for the 
fiber cross-section plane xy. The notations e r = e r (r), 
e- v = e ¥ ,(r), and e z = e z (r) stand for the cylindri¬ 
cal components of the profile function e^’+’+^r, p) of 
the forward counterclockwise polarized guided mode and 
are given in Refs. [HI, IHUHl ]. Equations m show that 
the x- and y-polarized guided modes have, in general, 
not only transverse but also longitudinal components. 
The local polarizations of these modes vary in the fiber 
cross-section plane, and are generally not strictly linear 
[HI [54 ]. It is interesting to note from Eqs. m that 
the signs of the longitudinal components fe z cos p and 
fe z sin p for the x- and y-polarized modes, respectively, 
depend on the propagation direction /. Thus, the differ¬ 


ence between the mode profile functions for the forward 
(/ = +) and backward (/ = —) guided fields is expressed 
by the change in sign of the longitudinal components. 
This change may affect the magnitude of the coupling 
between the atom and the field and, consequently, may 
lead to directional spontaneous emission and directional 

scattering [HI. [HI. l38l. IHI. IHl] . 

Since the radial direction of the atomic position in our 
study is parallel to the axis x, the polar angle for the 
atomic position is p = 0. Therefore, when we evaluate 
the mode functions m at the positions of the atoms, we 
find 

e (a,/x) | x axis = y/2(i\e r \i + /|e 2 |z), (2a) 

e ( “ /y) \ x axis = iV 2 \e v |y. (2b) 

In deriving the above equations we have used the prop¬ 
erties e r = *|e r |, e v = — |e^,j, and e 2 = \e z \ [HI IHl - l55j . 

Equation (I2al) indicates that, on the x axis, the x- 
polarized guided mode has two components, the trans¬ 
verse component iy/2\e r \ (aligned along the x axis) and 
the longitudinal component fy /2 |e 2 | (aligned along the 2 
axis). The difference in phase between these components 
is / 7t/2. It depends of the mode propagation direction /. 
The polarization vector of the ^-polarized guided mode 
on the x axis is 


i|e r |x + /|e 2 |z 
y/\e r \ 2 + | e 2 1 2 


It is clear that the local polarization of the rr-polarized 
guided mode on the x axis is elliptical in the zx plane. 
The ellipticity vector is given byijuxu*] = <7 e iiy, where 


Cell = / 


2|e r | |e 2 | 

M 2 + M 2 


( 4 ) 


is the ellipticity parameter. The circulation direction of 
the above elliptical polarization depends on the mode 
propagation direction /. The polarization vector ([3]) is 
a linear superposition of the circular polarization basis 
vectors er+ = (z + *x)/\/2 and cr~ = — (z — iic)/y /2 with 
the coefficients 

_ 1 |e r |+/|e 2 | 

Ua+ ~ v'Me.P + le^’ 

1 \e r \-f\e z \ U 

\/2 y/\e r \ 2 + |e 2 | 2 ' 


The corresponding weight factors |ucr±| 2 are given in 
terms of the polarization ellipticity parameter cr e ii as 
|m ct ±| 2 = (1 ± u e ii)/2. When |e 2 | ~ |e r |, we have 
(Tell — ± 1 , \u a ±\ 2 ~ 1, and |?Vt| 2 — 0 for / = ±. In 
this case, the local polarization of the ^-polarized guided 
mode at the position of the atom on the x axis is almost 
circular. As an example, we consider the case where the 
nanofiber radius is a = 250 nrn and the atom-to-surface 
distance is r — a = 200 nm. These parameters corre¬ 
spond to the Vienna atom trap experiment j35| . We find 









4 


the ratio |e z |/|e r | ~ 0.55, which leads to er c iL— ±0.84, 
\u a ±\ 2 ~ 0.92, and |w ct t| 2 — 0.08 for / = ± [33|, HU- 

Equation (I2bll indicates that, on the x axis, the y- 
polarized guided mode has a single component i\J2 je^l, 
which is aligned along the y axis. Thus, the local polar¬ 
ization of the y-polarized guided mode at the position of 
the atomic array is exactly linear along the y axis. This 
local polarization does not depends on the mode propa¬ 
gation direction /. 

B. Atomic levels and EIT schemes 

We assume that the atoms have a single upper level 
|e) of energy Tiui e and two lower levels | g) and \h) of en¬ 
ergies TiiOg and huih, respectively. The atoms are initially 
prepared in the lower level | g). The levels |e) and | g) are 
coupled by a weak probe field £ p of frequency ui p . The 
levels |e) and | h) are coupled by a strong control field £ c 
of frequency w c . The transition between the lower lev¬ 
els | g) and | h) is electric-dipole forbidden. We assume 
that the probe field £ p is guided by the nanofiber and 
propagates along the fiber axis ^ in the direction f p = ±, 
where + or — corresponds to the positive direction +z 
or the negative direction —z, respectively. The detuning 
of the probe field with respect to the atomic transition 
|e) -O- | < 7 ) is denoted by 

A = LOp - U] 0 , (6) 

where to 0 = ix eg = w e — uj g . The control field £ c is an 
external plane-wave field propagating perpendicularly to 
the fiber axis 2 . The two-photon (Raman) transition 
between the lower levels \g) and | h) may be off resonance, 
and the corresponding two-photon detuning is denoted 

by 

S = Up - UJ C - Ulh + U)g. (7) 

In our general analytical calculations, the single-photon 
detuning A and the two-photon detuning S can be dif¬ 
ferent from each other. However, in our numerical cal¬ 
culations, we will study only the case where A = <5, that 
is, the case where the control field is on exact resonance 
with the atomic transition |e) 0 | h). 

To be specific, we use the transitions between the Zee- 
man sublevels of the D 2 line of atomic cesium in our cal¬ 
culations. In order to specify the internal atomic states, 
we use the minor principal axis y as the quantization 
axis zq . The purpose of this choice is that it allows us to 
identify appropriate atomic states |e) and | g) which are 
coupled to each other by just one type of polarization 
of guided probe light. Indeed, as shown in the previous 
subsection, at an arbitrary position on the axis x, the 
local polarization of ^-polarized guided light is elliptical 
in the zx plane and the local polarization of y-polarized 
guided light is exactly linear along the y axis. Therefore, 
an atomic transition between two Zeeman sublevels spec¬ 
ified with respect to the quantization axis y can interact 


with either x- or y-polarized guided light but not with y- 
or :r-polarized guided light, respectively. 

We consider two schemes that are illustrated in parts 
(a) and (b) of Fig. [L] In both schemes, we use the excited- 
state sublevel \ 6 P 3 / 2 l F' = 4 ,M' = ±4) as the upper 
level |e). Furthermore, we use the ground-state sublevels 
16 S 1/2 ,F = 3 ,M = ±3) and \6S 1/2 ,F = 4 ,M = ±4) as 
the levels |<?) and | h), respectively, in scheme (a) of Fig. [I] 
and as the levels | h) and | g), respectively, in scheme (b) of 
Fig-Ill The effects of other Zeeman sublevels are removed 
by applying an external magnetic field. We emphasize 
that the atomic states in the above schemes are specified 
by using the minor principal axis y as the quantization 
axis. In addition, the quantum numbers of the lower 
states | < 7 ) and | h) are interchangeable between the two 
different schemes (see Fig. [T]). 

In scheme (a) of Fig. [TJ the guided probe field £ p is 
quasilinearly polarized along the x direction. Meanwhile, 
the control field £ c is a plane wave propagating along 
the axis x and linearly polarized along the axis y = zq 
in accordance with the 7 r type of the atomic transition 
|e) O | h). In scheme (b) of Fig. [0 the guided probe 
field £ p is quasilinearly polarized along the y direction. 
Meanwhile, the control field £ c is a plane wave prop¬ 
agating along the axis y and counterclockwise or clock¬ 
wise circularly polarized in accordance with the er + or <j_ 
type of the atomic transition |e) -H- | h). In what follows, 
schemes (a) and (b) of Fig. [T] are called the x- and y- 
polarization schemes, respectively. In both schemes, we 
neglect the reflection of the control field £ c from the fiber 
surface. Since the control field propagation direction and 
the atomic array axis are perpendicular and parallel, re¬ 
spectively, to the fiber axis z, the effect of the reflection of 
the control field £ c can be easily accounted for by mod¬ 
ifying the magnitude of £ c at the position of the atomic 
array. 

We introduce the notation M a for the magnetic quan¬ 
tum number of the atomic level a, where a = e, < 7 , h. 
In the analytical calculations, we use M g = ±3 and 
Mh = ±4 for the cc-polarization scheme, M g = ±4 and 
Mh = ±3 for the y-polarization scheme, and M e = ±4 
for both schemes. 

We note that, in both schemes (a) and (b) of Fig. [L] the 
probe transition |e) | g) is not coupled to the guided 

modes with the polarization that is orthogonal to the 
polarization £ p of the incident guided probe field. Here, 
we have introduced the notation £ = y for £ = x and 
£ = x for £ = y. For the control field £ Cl we can use 
a guided field instead of an external plane-wave field. 
Indeed, the control transition |e) | h) can be coupled by 

an additional guided field that is quasilinearly polarized 
along the y direction in scheme (a) of Fig. |T|or along the 
x direction in scheme (b) of the figure. 

Due to the interaction between the atoms and the 
guided probe field £ p , a guided reflected field £b with the 
frequency = u> p , the propagation direction /& = —f p , 
and the polarization £(, = may be generated. In order 
to describe the reflection, we need to include both prop- 


5 


agation directions f = + and / = — into the analysis. 
We introduce the notation £ f for the positive frequency 
component of the electric part of the field in the guided 
mode with the frequency w p , the polarization £ p , and the 
propagation direction / = ±. The electric field vector 
£/ = £f(r,if,z,t) is related to the photon flux ampli¬ 
tude Af = Af(z,t) via the formula 


£f=i ^P Af e^\ ( 8 ) 

V e 0 Vg 

Here, v g = (d/3/du})~ 1 \ ul=UJp is the group velocity of the 
guided field, where /? is the longitudinal propagation 
constant. The notation stands for the normal¬ 

ized profile function for the guided mode with the fre¬ 
quency oj, the propagation direction /, andMhe polar¬ 
ization £ [FI], [EH [HU ij56j]. According to jUj], the am¬ 
plitudes Af = A± of the photon fluxes of the guided 
fields £f = £± are governed, in the framework of the 
slowly varying envelope approximation, by the propaga¬ 
tion equations 

(J^ - Wp'j -A+ = n A G* + p eg , 

A- = - n A G*_p eg ■ (9) 

Here, /3 P = /3(w p ) is the longitudinal propagation con¬ 
stant for the forward and backward guided fields, 

n A = Y^Hz-Zj) ( 10 ) 

3 

is the one-dimensional atom-number density, G+ and G_ 
are the coupling coefficients for the guided fields £+ and 
£_, respectively, and p aa > = {a\p\a') with a,a' = e,g,h 
are the elements of the density matrix p the atom in the 
interaction picture. The coupling coefficient G/ = G± is 
defined as 


In the case of the ^-polarization scheme, the coupling 
coefficient Gf is given as 

Gf = MW + /(Me - Mg) |e,|). (12) 

Here, e r and e z are respectively the radial and axial 
components of the mode profile function e = e^ p,+,+ ^ 
of the forward counterclockwise quasicircularly polarized 
guided modes [5lL ;53j, [H, [Sbj . 

In the case of the y-polarization scheme, the coupling 
coefficient G/ is given as 

Gf=l \[^g de9KV (13) 

Here, e^ is the azimuthal component of the mode profile 
function e = e^ p,+ ’ + ) 0 f the forward counterclockwise 
quasicircularly polarized guided modes [5l|, [H, |H, @1 • 

Note that |G /| 2 = 7 eg^ P \ where 7 eg p ^ is the rate of 
spontaneous emission of the atomic transition |e) -f-t | g) 
into the guided mode with the propagation direction 
/ and the polarization £ p . Since the probe transition 
|e) O | g) is not coupled to the guided modes with the 
polarization £ p in the cases of schemes (a) and (b) of 

Fig. [Q we have the relation 7 iP = 7 e 9 ?p \ where "fiP is 
the rate of spontaneous emission of the atomic transition 
|e) O | g) into the guided modes with the propagation 
direction /. Hence, we obtain |G /| 2 = 7 iQ. 

The Rabi frequency caused by the guided field 5/ is 
given by 

flf = d eg • £f/h = 2iGjAf. (14) 

The Rabi frequency caused by the control field £ c is = 
d eh ■ £c/foi where d e h is the dipole matrix element for 
the atomic transition |e) o | h). We assume that |H±| -C 
|f2 c |. In the adiabatic approximation, the expression for 
p e g to first order in fl± is found to be 





( 11 ) 


Here, d eg is the dipole matrix element for the atomic 
transition |e) O | g). We emphasize that the photon flux 
amplitudes A±(z,t) are independent of x and y. 

It is clear from Eq. m that the atom-field coupling 
coefficient Gf depends on the local polarization of the 
guided probe field at the position of the atom. This coef¬ 
ficient also depends on the orientation and magnitude of 
the dipole matrix element d eg . We emphasize again that, 
in our study, the internal atomic states and the atomic 
transitions are specified by using the minor principal axis 
y as the quantization axis. Moreover, in order to obtain 
nonzero atom-field coupling coefficients, different transi¬ 
tions |e) 0 | g) of atomic cesium, with different dipole 
matrix elements d eff , are used in the different polariza¬ 
tion schemes (see Fig. [TJ) . 


where [lr:; 6 ] 


Peg = 7, -F 


T = 


5 + iThg 

W! 4 — (A + iT e g))8 + i£hg) 


(15) 


(16) 


Here, T eg is the decay rate of the atomic probe transition 
coherence p eg , and T h g is the decay rate of the lower-level 
coherence ph g - 

We insert Eq. m into Eqs. © and make use of ex¬ 
pression m- Then, we obtain 

A+ = i(K ++ A+ + Ah—*4_), 

+ *&) A - = -i(K-+A+ + K -4_), (17) 
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where 

Kff = n A G)G f ,T (18) 

for /, /' = +,—. In general, the one-dimensional atom- 
number density n A and consequently the coupling coeffi¬ 
cients I\ f ff are generalized functions of 2 [see Eq. ([TO ] . 

III. CONTINUOUS-MEDIUM 
APPROXIMATIONS 

In this section, we approximate the generalized- 
function representation m of the one-dimensional 
atom-number density n A by two different continuous- 
function representations and present the corresponding 
analytical and numerical results. In Sec. IIII Al we replace 
n A by a constant and solve the corresponding coupled¬ 
mode propagation equations. In Sec. IIII B1 we expand 
n A into a Fourier series and neglect the terms that do 
not correspond to the phase-matching condition in the 
coupled-nrode propagation equations. 

A. Homogeneous-medium approximation 

We consider the case where the lattice constant A is 
not close to any integer multiple of the in-fiber half¬ 
wavelength Aj?/2 = 7t//3 p of the probe field, that is, the 
atomic array is far off the Bragg resonance. In this case, 
the effect of the interference between the beams reflected 
from different atoms in the array is not significant and, 
therefore, we can neglect the discreteness and periodicity 
of the atomic array. This approximation means that we 
can use the one-dimensional atom-number distribution 

n A = 1/A, (19) 


condition is A-(L) = 0. In this case, the reflection and 
transmission coefficients are = Al_(0)/*4+(0) and 

= A + (L) / A + (Q), respectively. In the case where 
the incident probe field is A~(L), the boundary condition 
is -A+(0) = 0. In this case, the reflection and transmis¬ 
sion coefficients are R ^ = A+(L)/A~(L) and ^ = 
-4-(0)/„4-(A), respectively. With the help of the rela¬ 
tion K + _ = /i_ + , we can show that R ^ = R ^ = R A . 
The expressions for R A and T A are found to be 

R _ iK-+ sin (QL) _ 

A Q cos {QL) - i[/3 p + ( K + + AT_)/2] sin (QL) ’ 

T (f) = _ Qexp[if(K + - K-)L/2] _ 

A Qcos(QL) - i[j3 p + (K+ + A'_)/2] sin(QL) ’ 

where 

Q = \jvi + Pp(K+ +I<-)+ {K+ ~ S £. (23) 

We note that, in the case of the y-polarization scheme, 
we have K + = K_, which leads to = T A ^ =T A . 

Since \G^Gf>\ < 2T es , |< 1 /T eg , and A ~ ir/j3 p , we 
have \Kff | <C f3 p . Hence, we find Q ~ /3 p + (K + +K_)/2. 
With this approximation, Eqs. (l22l) reduce to 

Ra “ 2^nf7T7fr {exp[l(2 ' 3 >’ +K * + K - )L] ~ ll - 

T { A f) ~exp[i(p p + K f )L\. (24) 

Since | K _|_| <C /3 P , the reflectivity of the array is 

\R A \ 2 — 0. The transmission of the probe field is 
[T^l 2 ~ exp(—2 KfL), where 

K f = Im (K f ) (25) 


which is continuous and constant in the axial coordinate 
z. With the use of this approximation, expression m 
for the coefficients Kff reduces to 

Kff = Kff , ( 20 ) 


is the absorption coefficient for the probe field in the 
case where the reflection is negligible. The corresponding 
phase shift coefficient for the probe field is 

d f = Re(A7). (26) 


where the coefficients 


The optical depth per atom T>f = 2k f A is 


Kff = G}GfZ (21) 

are independent of z. Note that K H _ K _ = K ++ K _ 

In the case of the x-polarization scheme, we have A' + + 7 ^ 

A'__ and AT_|_= K _In the case of the y-polarization 

scheme, we have K ++ = AT__ = Ah| = AT It is 

convenient to introduce the notations K + = K ++ and 
K- = K _, 

We can easily solve Eqs. o with the constant coef¬ 
ficients given by Eq. We assume that z = 0 and 

z = L = {N — 1)A are the left- and right-edge posi¬ 
tions of the atomic medium, respectively. In the case 
where the incident probe field is ^4+(0), the boundary 


V f =2\G f \ 2 lm( 


5-NT 


hg 


|H c | 2 /4 - (A + zT e£ ,)(<5 + iT hg ) 


The phase shift per atom 0/ = Of A is 


(27) 


9/ |G/|2RC (|H c | 2 /4- (A + iTeg)(S + iT hg )} ‘ 

(28) 

We now describe time dependence of the guided probe 
field Af } where / = f p . In the case where the period of 
the atomic array is far from the Bragg resonance, the re¬ 
flection is, as shown analytically above and illustrated nu¬ 
merically in Figs. 0b) and 0b) below, negligible. Then, 
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in the frequency domain, the Fourier-transformed ampli¬ 
tude A f (z,u) = ( 2 - 7 r ) -1 / 2 A/(z, t')e i ( UJ ~ UJp ) t dt of the 

probe field is governed by the propagation equation 

dA f d Z ’ U) = + K f (u>)]A f {*,u)- ( 29 ) 

Here, we have introduced the notation Kf(u>) = Kff(u >), 
where Kff>{ui) = Kff>\ Up=ul is given by Eq. (|2H) with 
the substitution uj p = co. We neglect the dispersion of 

I 


the fiber-mode group velocity, that is, we take /3(ui) = 
/3 p + /3' p (uj — cjp), where u p is the central frequency of the 
input guided probe field. We expand Kf(ui) up to the 
second order of w — uj p as 

A/M = K f + K' f {u - co p ) + - Wp ) 2 , (30) 

where 


K _ \G f \ 2 _ 6 + iT hg _ 

f A |O c | 2 /4 - (A -(- *r eg )(<5 + iT hg ) ’ 

K , = \Gf\ 2 \n c \ 2 /A+(S + iT hg r 

f A [|^ c | 2 /4 - (A + + ir^)] 2 ’ 

X" — 2 [ A + ^ e 3 + 2(^ + ^r^ g )] |H e | 2 /4 + (<5 + iT hg ) 3 
V/_ A [|H c | 2 /4 - (A + iT eg )(5 + iT hg )) 3 


We substitute Eq. ( 15(1 into Eq. ( 1551 ) and perform the 
inverse Fourier transformation. Then, we obtain 


9A f (z,t) 

fdz 


= i{P p + K f )A f (z, t ) - (# + K' f ) dA f t ' t] 


; K f d 2 A f (z,t) 
2 dt 2 


(32) 


In general, Kf, K'f, and K'J are complex parame¬ 
ters. The imaginary and real parts of the parameter Kf, 
namely, the coefficients Kf = Im(A'/) and 9f = R e{I\f), 
are, as already discussed above, the absorption and phase 
shift coefficients, respectively, for the guided probe light 
field. We emphasize that the coefficients Kf, K'f, and 
K'j are the propagation characteristics for the photon 
flux amplitude Af(z,t). 

We analyze the case of exact one- and two-photon res¬ 
onances, that is, the case where A = 6 = 0. It is clear 
from the expression for I\f in Eqs. m that, when 
|fi c | 2 T hg T eg , the absorption coefficient Kf = Im(A'/) 
and the phase shift coefficient Of = Re(Aj) are small. 
These features are the signatures of EIT 1-Q. We note 
that the width of the corresponding transparency win¬ 
dow is given by Atamans = 1 /\JL Im(AT^). When the 
dephasing rate Th g is negligible, the expression for K" 
in Eqs. (1511) yields Im(AT') oc 7 eg^r es /A|r 2 c | 4 , leading to 

M 


AcJtrans CK 



(33) 


When the input probe pulse is long enough, the group 
velocity V g of the probe field is determined by the equa¬ 
tion 1/Vg^ = l/v g + Re(A'^). When the dephasing rate 


r^ g is negligible, the expression for K £ in Eqs. (1511) yields 
Re(ATy) oc "fiP/A\fl c \ 2 , leading to [l]-Q 


V g [f) 


1 Jep 
v g A \n c r 


(34) 


where the atom-number density ua in the framework 
of the homogeneous-medium approximation is given by 

Eq. (1151) . 

We calculate numerically the optical depth per atom 
Vf, the phase shift per atom 0/, the group-velocity 
reduction factor c/V g *\ the transmittivity |T^| 2 , and 
the reflectivity \Ra \ 2 as functions of the detuning A 
of the guided probe field. In the numerical calcula¬ 
tions presented in this paper, we use, as already stated, 
the outermost Zeeman sublevels of the hyperfine levels 
6 P 3 / 2 F' = 4, 6 Si/ 2 A = 3, and 6 Si/ 2 F = 4 of the D 2 line 
of atomic cesium, with the free-space wavelength Ao = 
852.35 nm. In the calculations for the i-polarization 
scheme, we use the levels |e) = | 6 P 3 / / 2 ,F' = 4 ,M' = 4), 

1 9) = |6 S 1/2 ,F = 3, M = 3), and | h) = |6 S 1/2 ,F = 
4, M = 4). In the calculations for the y-polarization 
scheme, we use the levels \e) = | 6 P 3 / 2 , F' = 4, M' = 4), 

I g) = |65i/2,P = 4 ,M = 4), and | h) = | 6 Si/ 2 ,F = 
3 ,M = 3). The rate r eg of decay of the atomic optical- 
transition coherence is calculated by using the results of 
Ref. (56|. The obtained value is F eg ~ 27r x 2.67 MHz. 
The corresponding value of the linewidth of the upper 
level |e) is 2r eg ~ 27 t x 5.34 MHz. This value is slightly 
larger than the literature value 70 ~ 27r x 5.2 MHz for the 
atomic natural linewidth [57], |58[ • The obtained increase 
of the atomic linewidth is caused by the presence of the 
nanofiber. The lower-level decoherence rate is assumed 
to be Th g = 27t x 50 kHz. This value is comparable to 
the experimental value of about 27rx 32 kHz, measured in 
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Probe field detuning N2n (MHz) 


FIG. 2. EIT of an ^-polarized guided probe field propa¬ 
gating along an array of cesium atoms with a cr-type probe 
transition. The optical depth per atom T>f (upper row), 
the phase shift per atom Of (middle row), and the group- 
velocity reduction factor c/Vg^ (lower row) are plotted as 
functions of the detuning A. The working levels are |e) = 
\ 6 P 3 / 2 ,F' = 4, M' = 4), | ff ) = | 6 S 1/2 ,F = 3, M = 3), and 
\h) = |65'i/2,T’ = 4, M = 4). The lower-level decoherence 
rate is = 2-7T x 50 kHz. The probe field propagates in the 
positive (left column) or negative (right column) direction of 
the fiber axis 2 . The fiber radius is a = 250 nm. The distance 
from the atoms to the fiber surface is r — a = 200 nm. The 
atom-number density is ua = 1/A with A = 498.13 nm. The 
control field is at exact resonance with the atomic transition 
|e) 44 |ft) and is a plane wave propagating along the x di¬ 
rection and polarized along the y direction. The intensity of 
the control field is I c = 1 mW/cm 2 (the corresponding Rabi 
frequency is fl c = O. 4670 ). 


the Vienna experiment [32|. The control field is at exact 
resonance with the atomic transition |e) 74 |ft). The fiber 
radius is a = 250 nm and the distance from the atoms to 
the fiber surface is r — a = 200 nm |35|. 

We plot in Figs. [2] and 0] the results of the calculations 
for scheme (a) of Fig. [Q where the guided probe field is 
x polarized. We show in Figs. [3] and [5] the results for 
scheme (b) of Fig. [L] where the guided probe field is y 
polarized. For the calculations of Figs. 0][5j we use the 
one-dimensional atom-number density ua = 1/A, where 
A = 498.13 nm. The chosen value of A is one-half of 
the in-fiber wavelength of the red-detuned standing-wave 
guided light field used in the Vienna atom trap experi¬ 


ment [35|. Such a value of the array period is far from 
the Bragg resonance (the Bragg resonant array period is 
A res = nA^/2, where n = 1,2,... and A p = 2 it//3 p ~ 
745.16 nm for the probe field with the free-space atomic 
resonance wavelength X p = Ao = 852.35 nm). 

Figures [2] and [3] show the familiar features of EIT: in 
the vicinity of the exact resonance, the optical depth T>f 
achieves a low minimum, the magnitude of the phase shift 
0 / is small but the slope of the phase shift is steep, and 

( f) 

consequently the reduction of the group velocity Vg is 


FIG. 3. EIT of a ff-polarized guided probe field propagating 
along an array of cesium atoms with a 7r-type probe transition. 
The optical depth per atom D (a), the phase shift per atom 0 
(b), and the group-velocity reduction factor c/V g (c) are plot¬ 
ted as functions of the detuning A. The working levels are 
|e> = I6-P3/2, F 1 = 4, M' = 4), |ff) = \ 6 S 1/2 ,F = 4, M = 4), 
and |ft) = |65 'i/2,F = 3 ,M = 3). The lower-level decoher¬ 
ence rate is F^g = 2 - 7 T x 50 kHz. The control field is at ex¬ 
act resonance with the atomic transition |e) <4 |ft) and is a 
plane wave propagating along the y direction with the coun¬ 
terclockwise circular polarization. The intensity of the control 
field is I c = 1 mW/cm 2 (the corresponding Rabi frequency is 
fl c = O.4370). Other parameters are as for Fig.0 The results 
do not depend on the propagation direction of the guided 
probe field. 


large [T]-0|. These features occur in a region of frequency 
which is narrower than the atomic natural lincwidth 70 
and is called the EIT window. 

Figures |D(a) and (5)(a) show that the transmittivity 
|T ^| 2 has a narrow peak at A = 0. We observe from 
Figs. Hb) and [51(b) that the reflectivity |AU | 2 is very 
small and is slightly asymmetric with respect to A. The 
asymmetry results from the fact that the relative phase 
/ 3 p L is not an integer multiple of n for the parameters 
used. Due to this fact, the magnitude of the reflectivity 
|I?a | 2 depends on the sign of the detuning A. 

When we compare the left and right columns of Fig. [2] 
where the probe field polarization is = x and the prop¬ 
agation direction is / = + for the left column and / = — 
for the right column, and compare these two columns 
with Fig. [3j where = y and / = ±, we see that the 
optical depth per atom Vf, the phase shift per atom Of, 
the group-velocity reduction factor c/Vg^ are different 
in the three cases. Similarly, when we compare the solid 
red and dashed blue curves of Fig.[4|a), where the probe 
field polarization is = x and the propagation direc¬ 
tion is f = + for the solid red curve and / = — for the 
dashed blue curve, and compare these two curves with 
the curve of Fig. [SJa), where (; p = y and / = ±, we ob¬ 
serve that the transmittivity |T ^| 2 is also different in 
the three cases. The differences between the results for 
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Probe field detuning A/2rt (MHz) 


FIG. 4. (Color online) Transmittivity \T^\ 2 (a) and re¬ 
flectivity |-Ra | 2 (b) of an x-polarized guided probe field as 
functions of the detuning A in the homogeneous-medium ap¬ 
proximation. The atom-number density is ua = 1/A with 
A = 498.13 nm. The medium length is L = (N — 1)A with 
N = 200. Other parameters are as for Fig. 0 



Probe field detuning AJ2n (MHz) 


FIG. 5. Transmittivity |Ta | 2 (a) and reflectivity |-Ra | 2 (b) 
of a y-polarized guided probe field as functions of the de¬ 
tuning A in the homogeneous-medium approximation. The 
atom-number density is tia = 1/A with A = 498.13 nm. The 
medium length is L = (N — 1)A with N = 200. Other param¬ 
eters are as for Fig. [3] 


the three cases arise from the fact that the atom-field cou¬ 
pling characterized by the coefficient Gf depends on the 
local polarization of the guided probe field at the position 
of the atom [see Eq. mi that is, on the mode polariza¬ 
tion £ p = x,y and on the propagation direction / in the 
case of = x [see Eqs. CE2D and mi Another reason 
is that different atomic transitions |e) |g), which have 

different complex dipole matrix elements d eg , are used in 
the x- and y-polarization schemes (see Fig. [[]). 

Thus, in the case of the x-polarization scheme, the op¬ 
tical depth Vf, the phase shift 0 /, the group velocity 


Vg^\ and the transmittivity \T^\ 2 substantially depend 
on the propagation direction /. The directionality of 
transmission of guided light through the array of atoms 
is a consequence of the existence of a longitudinal com¬ 
ponent of the guided light field as well as the ellipticity of 
both the field polarization and the atomic dipole vector 
f5ll . fbbf . We note that directional spontaneous emission 
into an optical nanofiber has been recently demonstrated 
experimentally for trapped atoms [33| and for nanopar¬ 
ticles [38|. An optical diode based on the chirality of 
guided photons has also been reported [34] . 

We observe from Figs. UKa), 0[b) , and0](a) for the x- 
polarization scheme and Figs. [5](a) and [SJa) for the y- 
polarization scheme that, at zero detuning, the absorp¬ 
tion is not completely suppressed and, hence, the induced 
transparency is not 100%. This is a consequence of the 
fact that a nonzero lower-level decoherence rate Th g is 
included in our calculations. Without this decoherence, 
the induced transparency is 100% at zero detuning. Com¬ 
parison between Figs. [21(a), 0(b), and 0(a) and between 
Figs. SK a ) andO a ) shows that, for a given lower-level de¬ 
coherence, the residual absorption depends on the probe 
field polarization £ p and on the propagation direction / 
in the case of £ p = x. The reason is twofold: (1) the 
atom-field coupling depends on the local polarization of 
the field and ( 2 ) the different atomic transitions |e) F 4 | g) 
are used in the different polarization schemes. 


B. Phase-matching approximation 

We now take into account the periodicity of the 
atomic array but still consider the array as a contin¬ 
uous one-dimensional medium. For this purpose, we 
rewrite expression mi for the one-dimensional atom- 
number distribution ua in the Fourier series form ua = 
(i/A) E^oc e 2mriz / A for 0 < z < L = {N- 1)A. Then, 
Eq. 1151) for the coupling coefficients Kff becomes 

OO 

Kff = Kff E e2ni01 * tZ • ( 35 ) 

n=— oo 

Here, the parameter /?i at = n/A is the wave number that 
characterizes the periodicity of the atom distribution in 
the array. The coupling coefficients Kfy are independent 
of z and are given by Eq. m- 

We assume that there exists a positive integer num¬ 
ber n = 1,2,... such that the propagation constant /3 P 
for the guided probe field is close to n/3i at . This spe¬ 
cific integer number n indicates the dominant role of the 
corresponding harmonic in Eq. (1351) . We note that the 
equality /3 p = nf 3i a t means /3 P A = nir, which is the geo¬ 
metric condition for the nth-order Bragg resonance. 

We introduce the transformation A+ = a+e m/3lat2 
and A- = a_e -m ^ lat2 . When we insert Eq. (1551) into 
Eqs. da and keep only the phase-matching terms, we 
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obtain 

( d \ 

f g - i5/3 p J a + = il\ ++ a + + iK + _a-, 

—f i8/ 3p^j a- = —iK — |_a+ — z/v__a_. (36) 

Here, the factor S/3 P = (3 P — n/ 3i at characterizes the phase 
mismatch between the guided probe field and the nth 
harmonic of the atomic distribution. In Eqs. (1551) . the 

terms iK ++ a+ and — iK _o_ describe the self coupling 

of the modes and originate from the zerotli-order har¬ 
monic in the Fourier expansion of n a, while the terms 

iK. (_a_ and — iK _|_a+ describe the coupling between 

the forward and backward fields and originate from the 
nth-order harmonic in the Fourier expansion of ua ■ Sim¬ 
ilar to the results of the previous subsection, the reflec¬ 
tion and transmission coefficients in the phase-matching 
approximation are found to be 

zA'_ + sin(17L) 

A U cos (UL) — i[S(3 p + (K + + A'_)/2] sin((7L) ’ 

T U) = U ex P [if( K + ~ K-)L/2] exp(w/3i a tL) 

A U cos(UL) — i[5/3 p + (K + + K_)/2] sin(I7 L)' ’ 

where 

u=\jmr + 6P P (K+ + K-)+ (A+ ~ A ~ )2 . (38) 



Probe field detuning A/2 ti (MHz) 


FIG. 6. (Color online) Transmittivity IT^I 2 (a) and reflec¬ 
tivity |-Ra | 2 (b) of an s-polarized guided probe field as func¬ 
tions of the detuning A in the phase-matching approximation. 
The array period is A = 498.13 nm and the array length is 
L — (N — 1)A, where N = 200. Other parameters are as for 
Fig. [U 

The condition for the validity of the above results is 
\5(3 p \ <C /3 P . When the phase-mismatch factor S/3 P is 
sufficiently large that \Kfp\ <C |5/3 P |, Eqs. (1371) reduce 
to 

Ra ~ 2^+ A g( + A_ teXp|,(2 ^ + IU + K ~ )L] - ll ' 
ry>~exp[i((3, + A7)L]. (39) 



FIG. 7. Transmittivity |Ta| 2 (a) and reflectivity \Ra\ 2 (b) of 
a j/-polarized guided probe field as functions of the detuning 
A in the phase-matching approximation. The array period is 
A = 498.13 nm and the array length L = (N — 1)A, where 
N = 200. Other parameters are as for Fig. [3] 

Under the condition 1 1\ _|_| <C the reflectivity 

of the array is |-Ra| 2 — 0, while the transmittivity is 
|Ty } | 2 ~ exp(—2 k/L), with the absorption coefficient 
k f being given by Eq. (1751) . 

We plot in Figs. [G] and [7] the tuning dependences of 
the transmittivity \T^\ 2 and the reflectivity |IiU| 2 , cal¬ 
culated from Eqs. (T37l) for the parameters of Figs. 0] 
and[5j We observe that Figs. [51(a) and Ufa) are almost 
identical to Figs. Oja) and Eta), respectively. However, 
the peak values of the curves in Figs. [51(b) and [3(b) are 
one order of magnitude larger than that in Figs. mb) 
and 0Kb), respectively. The discrepancy is due to the 
fact that the periodicity of the atomic array is neglected 
in the homogeneous-medium approximation but is taken 
into account in the phase-matching approximation. 


IV. TRANSFER MATRIX FORMALISM FOR A 
DISCRETE ARRAY 

We now take into account not only the periodicity of 
the atomic position distribution ua = ~~ z j ) but 

also the discrete nature of this distribution. For this pur¬ 
pose, we use the transfer matrix formalism. 

A. Transfer matrix and input-output relation 

We first consider atom j with the axial coordinate z :j . 
We introduce the notations = lim^o.,- (zj ± e) for 
the right- and left-hand-side limiting points. The fields 
A+ {zj ) and A _ (Zj ~) with the propagation directions / = 
+ and / = —, respectively, at the limiting points z~ and 
Zj, respectively, can be interpreted as incoming fields 
with respect to the atom. Meanwhile, the fields A+(zj~) 
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and A- (Zj ) with the propagation directions / = + and 
/ = —, respectively, at the limiting points and z~, 
respectively, can be considered as outgoing fields with 
respect to the atom. According to the causal principle, 
the atom interacts with the incoming fields A±(z^f) but 

I 

not with the outgoing fields A±(Zj ). 

We introduce the notations Al = A{z^) and Ar = 
Aizjf), where A is the vector consisting of the com¬ 
ponents Af. When we integrate Eqs. (1171) over an in¬ 
finitely small interval dz around the position Zj of the 
atom and follow the procedures of Ref. [HlJ, we find the 
input-output relation 


Ar = WAl, (40) 

where 

W = T Ar_1 M (41) 


is the transfer matrix for the whole atomic array, with 

T = MF (42) 


being the transfer matrix for a single spatial period of 
the array. The matrix M in Eqs. m and H 2 l) is the 
transfer matrix for a single atom. The matrix elements 
of this matrix are given as 


— 1 — S++ + 


S+-S-+ 

1 + S— ' 


M _= 


1 


Mj__ = - 


1 + S— 
S+- 


M_ + = - 


1 + 5 __ 
5-+ 

1 + S—' 


(43) 


where 


S fr = -ifAKff. = —ifG*fG f 'F. (44) 

Note that we have S.\ = — S-+ and S ++ S = 

S _|_ S _The matrix F in Eq. flUD is the atom-free 

guided-fielcl propagator. The elements of this matrix are 
given as 


Fff = e i ^ pA 6 f fr. (45) 


It is easy to diagonalize the single-period transfer ma¬ 
trix T. Using the result of this diagonalization, we find 
the following expressions for the elements of the transfer 
matrix W for the array: 


W ++ = 


z n/2 f M++ sin(lVC) _ .,i fl „ A sin[(iV- 1)C] ) 


l Vz sin C 


sin C / 


w __ = z n/2 f M — sin(A^C) _ . ifl „A sin[(JV- 1)C] 1 
1 a/Z sin £ sin £ / ’ 

W_ + = -W + _ = z n/2 M-+sHNQ . 


s[Z sin £ 


(46) 


Here, we have introduced the notations 

1-S++ _ l + i\G+\ 2 T 


Z = det(M) = 


1 + 5_ l + i\ G-\ 2 T 


and 


where 


C = i ln(D ± i\J 1 — D 2 ), 


D = —^=(Af ++ e i/3pA + M _ e~ i P pA ). 

2v Z, 


(47) 


(48) 


(49) 


Note that £ is, in general, a complex parameter and sat¬ 
isfies the relations cos £ = D and sin £ = + \/l — D 2 . 

In the case of the ^-polarization scheme, the elements 
of the transfer matrix M for a single atom are found to 
be 


M ++ = 1 - 
M__ = 


S r 


2 S r 


1- S r ~ S z 

1 


2 S r 


1 - S r - S z 

M+_ = -M_+ = 


2 Srz 
S r -S z 


1 — S r ~ S z + 2 S rz 

Here, we have introduced the notations 


(50) 


S r = -i^ZL\e r \ 2 T, 

2e 0 hv g 

S = _.Up 4 g_ 2 T 


Srz — 


2e 0 hv g 

°pd e g 

'2eohv g 


(M e - M g )\er\\e z \T. 


(51) 


We find Z 7 +I in the case of the x-polarization scheme. 

In the case of the y-polarization scheme, the elements 
of the transfer matrix M for a single atom are found to 
be 

M++ = 1 “ 2S{p 


M _= 


1-5. 


1 -S v ’ 


M+- = —M-+ = — 


1 — 5. 


Here, we have introduced the notations 

^Kf 

e 0 hvg 


(52) 


(53) 


We find Z = 1 in the case of the y-polarization scheme. 


B. Reflection and transmission of probe light 

For the guided probe field propagating in the direc¬ 
tion / = +, the reflection and transmission coefficients 
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of a single atom are given by R ^ = —M ( ,/M and 

T = Z/M__, respectively. The explicit expressions 
for these coefficients are 

R {+) = S-+, 

T (+) = 1 - S++. (54) 

The corresponding reflection and transmission coef¬ 
ficients of the atomic array are given by R^ = 

—W-+/W and = Z N /W , respectively. 

For the guided probe field propagating in the direction 
/ = —, the reflection and transmission coefficients of a 

single atom are given by R^> = M _j_ /M _and = 

respectively. The explicit expressions for these 
coefficients are 

R (_) = - 5 +_, 

= 1 + S _ (55) 


The corresponding reflection and transmission coef¬ 
ficients of the atomic array are given by R^ N ^ = 

W+-/W _and = 1 /W__, respectively. 

It follows from the properties 5_|_= —S _and 

W .1 = —W that the single-atom reflection coefficient 

and the linear-array reflection coefficient R ^ do 
not depend on the field propagation direction /, that 
is, R = R= R and R^ = R^ = Rn- However, 
the corresponding transmission coefficients and 
may depend on the propagation direction /. We find 


Rn = 


'T'(f) _ 

1 N ~ 


Rsm(N() 


sm(N() - v / 2 : T(-)e i ^ A sin[(7V - 1)C] ’ 
Z(i+/ Ar )/ 2 T(-)sinC 
sin(iV£) — \/~Z T(~'>e l P pA sin[(7V — 1)£] 


(56) 


( f ) 

We can show that Rn and Tjj satisfy the recurrence 
formulas 


Rn+i 


Rn + 


T^T^Re 2 ^ 

1 - R N Re 2i ^ A ’ 


T Cf) T (/) e ^ A 
W+1 1 - R N Re 2i P pA ' 


(57) 


which are in agreement with ray optics. It is clear that 
the reflection and transmission coefficients Rn and 
for the atomic array depend on the polarization of the 
guided probe light field. We note that, in the vicinity of 
a Bragg resonance, we can reduce Eqs. (15(1 to Eqs. (1371) . 
which were obtained in the phase-matching approxima¬ 
tion. 

The phases of the transmission and reflection coef¬ 
ficients are given by 4 >^' 1 = Im(lnTjy^) and = 
Im(lnJ?jv), respectively. The group delays of the trans¬ 
mitted and reflected fields are given by and 

T R = respectively, where = d&/doj p . Thus, from 
the frequency dependences of the transmission and re¬ 
flection coefficients we can calculate the group delays of 
the transmitted and reflected fields. 


It follows from the expression for in Eqs. (l56l) that 


T 


(+) 


N 


T, 




= z N = 


N 


1 + i\G+\ 2 T 
l + i\G-\ 2 T 


N 


(58) 


In the case of the x-polarization scheme, we have |G+| ^ 
|G_| and, hence, Z / 1. In this case, the single-atom 
and linear-array transmission coefficients T(T) and 
depend on the field propagation direction /. In the case 
of the y-polarization scheme, we have G+ = G_ and, 
hence, Z = 1. In this case, the single-atom transmission 
coefficient T^\ the linear-array transmission coefficient 
Tjy, and the group delay are independent of the 
propagation direction /, that is, we have T = T^~> = 
T, T^ +) = T^ -) = T n , and r^ +) = r!f ] = t t . 

The dependence of the single-atom transmission coef¬ 
ficient on the propagation direction / = ± in the 
case of x-polarized guided light is a consequence of the 
difference between the coupling coefficients G + and G_ 
for the different propagation directions + and —, respec¬ 
tively. The independence of the single-atom reflection 
coefficient R from the field propagation direction / is 
a consequence of the fact that R is proportional to the 
product of both coupling coefficients G + and G_. The 
dependence of the linear-array transmission coefficient 
Tx on the propagation direction / in the case of de¬ 
polarized guided light is a consequence of the difference 
between and for each atom in the array. The 
independence of the linear-array reflection coefficient Rn 
from the field propagation direction / is a consequence 
of the fact that all the atoms in the array considered here 
have the same reflection coefficient R. We note that the 
above properties of reflection and transmission of guided 
light in the atomic array are different from that of light 
in a conventional Fabry-Perot cavity formed by two mir¬ 
rors, with the reflection coefficients R[ + ^ = R[~^ = R\ 
and I? 2 + ^ = R-i 1 — R‘i an d the transmission coefficients 
T[ +) = T 1 (_) = Ti and T 2 (+) = T 2 (_) = T 2 . The transmis¬ 
sion coefficient T cav of such a cavity does not depend on 
the propagation direction. However, in the case where 
R\ ^ R 2 or Ti ^ T 2 , the cavity reflection coefficient R^H 
may depend on the propagation direction /. 

We calculate numerically the transmittivity \T^\ 2 , 
the reflectivity |i?Ar| 2 , the transmitted-field group delay 
Trp ', and the reflected-field group delay tr as functions 
of the detuning A of the guided probe field. We plot in 
Figs. [ 8 ] and [9] the results for the x-polarization scheme, 
and in Figs. [TO] and [Tl] the results for the y-polarization 
scheme. In the calculations for these figures, we used 
the value A = 498.13 nm, which corresponds to the ar¬ 
ray period in the situation of the atom trap experiment 
[35| ]. This value of the array period is, as already noted 
in the previous section, far from the Bragg resonance. 
We observe that Figs. OH a) and [TO^a) are almost iden¬ 
tical to Figs. [D(a) and Ola), respectively. This means 
that the homogeneous-medium approximation is a very 
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Probe field detuning A/271 (MHz) 


FIG. 8. (Color online) Transmittivity \T^\ 2 (a) and reflec¬ 
tivity |i?jv| 2 (b) of an a;-polarized guided probe field as func¬ 
tions of the detuning A. The array period is A = 498.13 nm, 
which is far from the Bragg resonance. The atom number is 
N = 200. Other parameters are as for Fig. [2] 



Probe field detuning A/2n (MHz) 

FIG. 9. (Color online) Group delays r ^ (a), rfjT^ (b), and 
tr (c) of an x-polarized guided probe field as functions of the 
detuning A. The array period is A = 498.13 nm, which is 
far from the Bragg resonance. The atom number is N = 200. 
Other parameters are as for Fig. [2] The dotted red line is for 
the zero group delay and is a guide to the eye. 


good approximation for the calculations of the transmit¬ 
tivity in the case where the array period is far from the 
Bragg resonance. Figures [SJb) andflUlb) show that the 
discrete-array reflectivity |i?Ar| 2 , like the homogeneous- 
medium reflectivity \Ra\ 2 hr Figs. 0b) and0b), is very 
small and is slightly asymmetric with respect to A. How¬ 
ever, the maximum magnitude of discrete-array reflec¬ 
tivity \Rn\ 2 is about twenty times larger than that of 
the homogeneous-medium reflectivity |i?^| 2 - Thus, the 
homogeneous-medium approximation is not reliable in 
evaluating the reflectivity. The discrepancy is due to the 
fact that the reflected field results from the interference 
between the fields reflected from different atoms in the 



Probe field detuning A/2jt (MHz) 


FIG. 10. Transmittivity |T)v| 2 (a) and reflectivity |-Rjv| 2 
(b) of a y-polarized guided probe field as functions of the 
detuning A. The array period is A = 498.13 nm, which is 
far from the Bragg resonance. The atom number is N = 200. 
Other parameters are as for Fig. [3] 



Probe field detuning A/2n (MHz) 


FIG. 11. (Color online) Group delays tt (a) and tr (b) of 
a 7/-polarized guided probe field as functions of the detuning 
A. The array period is A = 498.13 nm, which is far from 
the Bragg resonance. The atom number is N = 200. Other 
parameters are as for Fig. [3] The dotted red line is for the 
zero group delay and is a guide to the eye. 

array. The periodicity of the array greatly affects the in¬ 
terference but is neglected in the homogeneous-medium 
approximation. Comparison between Figs. 0b), 0b), 
and0b) and between Figs. 0b), 0b), and [TUI b) shows 
that the reflectivity calculated in the phase-matching ap¬ 
proximation is closer to the rigorous result |i?/v| 2 than 
that calculated in the homogeneous-medium approxima¬ 
tion. 

Figures0a),0b), and fTTT a) show that the group delay 
( f) 

Trp ' of the transmitted field has a positive peak at A = 0. 
This result means that the transmitted field is slow light 
in the vicinity of the atomic resonance. Close inspection 
shows that the results for the group delay r^ plotted 
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in Figs. Ha), [lb), and II If a) are in full agreement with 
the results for the group-velocity reduction factor cjV^ ] 
plotted in Figs. He), 0 f), and He), respectively. Ac¬ 
cording to Figs. H c ) andQjJb), the group delay tr of 
the reflected field has a negative-valued dip at A = 0. 
This result means that the reflected field is fast light in 
the vicinity of the atomic resonance Ili0. Com¬ 
parison between the solid red and dashed blue curves of 
Fig. Ha) and between Figs.H a ) andHb) shows that the 
transmittivity \T^\ 2 and the corresponding group delay 
Tjf) depend on the propagation direction / in the case of 
the ^-polarization scheme. This dependence is a conse¬ 
quence of the directional dependence of the coupling be¬ 
tween the z-polarized guided probe field and the atomic 
transition |e) | g), which is of the a type with respect 

to the y axis in the case considered [see Eq. (TH21) ] . We ob¬ 
serve from Figs.lQlandfTTlthat the dip of the reflected-field 
group delay tr has a narrower width and a larger magni¬ 
tude than the peak of the transmitted-field group delay 
Tj, . The difference is due to the fact that the transmit¬ 
ted field is determined by the interference between the 
incident field and the single-atom scattered fields while 
the reflected field is determined by the interference be¬ 
tween the fields reflected from different atoms. 



Time t (ps) 


FIG. 12. (Color online) Time dependences of the normal¬ 
ized intensities |Aout,ref/Ao| 2 of the transmitted and reflected 
pulses in the case of the ^-polarization scheme. The input 
pulse is of the Gaussian form, with the central frequency 
CJ P = too, the pulse length to = 2 fis, and the peak time 
t = 0. The number of atoms is N = 200. The array period 
is A = 498.13 nm. Other parameters are as for Fig. [2] The 
dotted red line is for the input pulse peak time t = 0 and is a 
guide to the eye. 


C. Time evolution of guided probe field pulses 

We now consider the time evolution of guided field 
pulses propagating along the atomic array. We as¬ 
sume that the guided input pulse is incident onto 
the atomic array in the / direction. The ampli¬ 
tudes -4i„(f) and ii n (u) of the incident field in the 
time and frequency domains, respectively, are related 
to each other by the Fourier transformation Ai n (w) = 
(27 t) -1 / 2 A ln (t)e lujt dt. The time-dependent trans¬ 
mitted field is 

A out (t) = — T^\u)AMe~ iult du:, (59) 

V 2-7T J — oo 

where T^\u>) = \ Up =ui- The time-dependent re¬ 
flected field is 

1 r°° 

Arei(t) = —== / Rn( u})Ai n (to)e iUt duo , (60) 

V 27T J — oo 

where R n (u>) = Rn\u> p =u>- 

We plot in Figs. OH and OH the time dependences of 
the normalized intensities |-A ou t,ref/dlo | 2 of the trans¬ 
mitted/reflected pulses in the cases of the x- and y- 
polarization schemes, respectively. In the numerical cal¬ 
culations for these figures, we used an input pulse of the 
Gaussian form Ai n (t) = Aoe _t / T oe _lajj,t , where r 0 is the 
initial characteristic pulse length. The Fourier transform 
of the input pulse is Ain(w) = (Aqt 0 / v / 2)e _T o ^ _ajp ^ 2 / 4 . 
We observe that, when the probe pulse is incident onto 



Time t (ps) 


FIG. 13. (Color online) Time dependences of the normal¬ 
ized intensities |A 0 ut,ref/Ao| 2 of the transmitted and reflected 
pulses in the case of the y-polarization scheme. The input 
pulse is of the Gaussian form, with the central frequency 
LO p = too, the pulse length ro = 2 fis, and the peak time 
t = 0. The number of atoms is N = 200. The array period 
is A = 498.13 nm. Other parameters are as for Fig. [3] The 
dotted red line is for the input pulse peak time t = 0 and is a 
guide to the eye. 


the array in the direction / = + in the case of the x- 
polarization scheme [see Fig. [T2f a)] and when the probe 
pulse is incident onto the array in an arbitrary direc¬ 
tion / = ± in the case of the y-polarization scheme [see 
Fig. QH a )], the transmitted pulse is weakened and slowed 
down significantly by the atomic array. For N = 200 
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atoms in the array with the length L = (N — 1)A ~ 100 
/rm, we obtain the group delays of about 1.17 p,s in the 
case of Fig. H2Ta) and 0.53 ps in the case of Fig. fldlab 
These group delays correspond to the group-velocity re¬ 
duction factors c/Vg ~ 3.5 x 10 6 and ~ 1.6 x 10 6 in the 
cases of Figs. fl 2 l a) and IT3l a) . respectively. Figure fT2lbl 
shows that, when the probe pulse is x polarized and inci¬ 
dent onto the array in the direction / = —, the intensity 
reduction and the group delay of the transmitted pulse 
are not significant. According to Figs. fl 2 Tcl and fl^l bl . 
the reflected pulse has a negative group delay (fast light) 
snun®. However, the intensity of this fast light is 
negligible. 


D. Bragg resonance 

We examine the transmission and reflection of the 
guided probe field under the EIT condition in the spe¬ 
cific case where the geometric Bragg resonance condition 
is satisfied, that is, j3 p A = nir , with n = 1,2,... being the 
order of the Bragg resonance. This condition involves the 
frequency lo p of the guided light field. Therefore, when 
we vary ui p in an interval A uj p around a Bragg resonance 
frequency, the Bragg resonance condition j3 p A = nn will 
be broken. However, if the frequency variation interval 
A uj p is small as compared to the optical frequency uj p , 
the effect of the deviation from the Bragg resonance can 
be neglected. 


1. Dependence on the probe field detuning 

We plot in Figs. [TTlfTTI the transmittivity |T^| 2 , the 

reflectivity |i?jv| 2 , the transmitted-field group delay t^\ 
and the reflected-held group delay tr as functions of 
the detuning A of the guided probe held in the case 
where the array period is A = 745.16 nm. This value 
of the array period satisfies the second-order Bragg res¬ 
onance condition A = A res = ttAf/2, where n = 2 and 
Xp = 27 t//3 p ~ 745.16 nm for the probe held with the 
atomic resonance wavelength X p = Ao = 852.35 nm. We 
avoid the hrst-order Bragg resonance with the aim of 
minimizing the effects of the direct dipole-dipole inter¬ 
action between the atoms. It is worth mentioning here 
that, in the framework of our treatment, the magnitudes 
of the reflectivity |i?jv | 2 and transmittivity |T/v | 2 of the 
atomic array for the guided helds do not depend on the 
order of the Bragg resonance. 

Figures [14] and [15] show the numerical results for the x- 
polarization scheme. Comparison between Fig. 1141 where 
the array period A is in the Bragg resonance, and Fig. [HI 
where A is far from the Bragg resonance, shows that the 
positive-direction (/ = +) transmittivity |T ^ f '*| 2 is al¬ 
most the same in the two cases [see the solid red lines in 
Figs. [bHal and[51[a)]. Meanwhile, the negative-direction 
(/ = —) transmittivity |T^A | 2 in the Bragg-resonance 
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Probe field detuning A/2 ji (MHz) 

FIG. 14. (Color online) Transmittivity |T^| 2 (a) and re¬ 
flectivity |_Rjv| 2 (b) of the atomic array in the axpolarization 
scheme as functions of the detuning A. The period of the ar¬ 
ray is A = 745.16 nm, which satisfies the second-order Bragg 
resonance condition. The number of atoms in the array is 
N = 200. Other parameters are as for Fig. [2] 



Probe field detuning A/2n (MHz) 

FIG. 15. (Color online) Group delays (a), (b), and 

tr (c) in the axpolarization scheme as functions of the detun¬ 
ing A. The period of the array is A = 745.16 nm, which sat¬ 
isfies the second-order Bragg resonance condition. The atom 
number is N = 200. Other parameters are as for Fig. [2] The 
dotted red line is for the zero group delay and is a guide to 
the eye. 


case is higher than that in the far-off-Bragg-resonance 
case, except for a very narrow frequency region around 
A = 0, where it is the same [see the dashed blue lines in 
Figs. [14K a) and[HKa)]. The reflectivity \Rn\ 2 in the Bragg- 
resonance case [see Fig. [T4l b)] is dramatically higher than 
that in the far-off-Bragg-resonance case [see Fig. [5Kb)]. 
The dependence of |l?jv | 2 on A is almost symmetric in 
the case of the Bragg resonance [see Fig. ICT bl]. Com¬ 
parison between Figs. EDA) and EKa) shows that the 
positive-direction transmitted-field group delay r ^ is 
almost the same in the two cases. Comparison between 
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flOl' b') shows that, due to the Bragg resonance condition, 
the reflectivity \Rn\ 2 is dramatically increased. 


2. Dependence on the atom number 

We plot in Figs. [18] and [T9] the dependences of 
the transmittivity |T^| 2 , the reflectivity |-Rat| 2 , the 
transmitted-field group delay , and the reflected-field 

group delay tr on the atom number N for the array 
period A = 745.16 nm, which satisfies the second-order 
Bragg resonance condition. We assume the atomic reso¬ 
nance A = 0 in these calculations. 


FIG. 16. Transmittivity \Tn\ 2 (a) and reflectivity |Ajv | 2 (b) 
of the atomic array in the y-polarization scheme as functions 
of the detuning A. The period of the array is A = 745.16 nm, 
which satisfies the second-order Bragg resonance condition. 
The atom number is N = 200. Other parameters are as for 
Fig. E] 



Probe field detuning A/2 ji (MHz) 



Atom number N 


FIG. 18. Dependences of the transmittivities |T^| 2 and 
\T^ ^| 2 , the reflectivity |-Rzv| 2 , and the group delays r^\ 
t^~\ and tr on the atom number N in the case of the x- 
polarization scheme. The period of the array is A = 745.16 
nm, which satisfies the second-order Bragg resonance con¬ 
dition. The detuning of the probe field is A = 0. Other 
parameters are as for Fig. [2] 


FIG. 17. (Color online) Group delays tt (a) and tr (b) in the 
2 /-polarization scheme as functions of the detuning A. The 
period of the array is A = 745.16 nm, which satisfies the 
second-order Bragg resonance condition. The atom number 
is N = 200. Other parameters are as for Fig. [3] The dotted 
red line is for the zero group delay and is a guide to the eye. 


Figs, m b) and [9[b) and between Figs. ITSl c) and [9[c) 
shows that the Bragg resonance modifies the values of the 
negative-direction transmitted-field group delay and 
the reflected-field group delay tr. However, the changes 
are not dramatic. 

Figures [16] and |TT] show the results of the numerical 
calculations for the y-polarization scheme. Comparison 
between Figs. fTHT a) audfTUTa) and between Figs. [IT] and 
HH shows that the effects of the Bragg resonance on the 
transmittivity |Tjv| , the transmitted-field group delay 
tt, and the reflected-field group delay tr are noticeable 
but not dramatic. Comparison between Figs. 116( b) and 


Figure [18] shows the numerical results for the x- 
polarization scheme. We observe from the left column 
of Fig. [18] that, when N increases, the positive-direction 
transmittivity |T^| 2 decreases to zero, the negative- 
direction transmittivity | Ty •*| 2 decreases to a nonzero 

limiting value iTicT^ 2 ~ 0.84 [see Eq. (I62bl) in the next 
part], and the reflectivity |l?Ar| 2 increases to a limit¬ 
ing value |i?oo| 2 — 0.08, which is smaller than unity 
[see Eq. m in the next part]. Figures [THKb) and 
IT8fd) show that, when N increases, the positive-direction 
transmitted-field group delay r^ + ' increases almost lin¬ 
early and hence can achieve an arbitrarily large value, 
while the negative-direction transmitted-field group de¬ 
lay ■* increases to about 0.1 /rs and then decreases in 
the range N < 4000 of the figure. Additional calcula¬ 
tions for N > 4000 show that ■* decreases to about 
24 ps at N ~ 7400 and then starts to slowly increase. 
We observe from Fig. 11514 ) that the reflected-field group 

































17 



Atom number N 


FIG. 19. Dependences of the transmittivity |T/v| 2 , the re¬ 
flectivity |i?jv| 2 , and the group delays tt and tr on the atom 
number N in the case of the y-polarization scheme. The pe¬ 
riod of the array is A = 745.16 nm, which satisfies the second- 
order Bragg resonance condition. The detuning of the probe 
field is A = 0. Other parameters are as for Fig. [3] 


delay tr is negative and its absolute value \tr\ decreases 
with increasing N in the range N < 4000 of the figure. 
Additional calculations for N > 4000 show that tr be¬ 
comes positive at N ~ 8440 and then tends to increase 
very slowly. 

Figure [!!)] shows the numerical results for the y- 
polarization scheme. We observe from the figure that, 
when N increases, the reflectivity \Rn\ 2 increases and 
approaches unity [see Eq. (165al) in the next part], the 
transmitted-field group delay tr approaches a limiting 
value of about 3 //s, and the reflected-field group delay 
tr is negative and approaches zero. 


|e r | > |e z | > 0, we have |i?oo| 2 < 1, that is, the limiting 
value |i?oo| 2 of the reflectivity for the ^-polarized guided 
fields is strictly smaller than unity. It is interesting to 
note that Eq. (RJTl) coincides with the result of Ref. [551 
for the case where the initial state of the atoms is an 
incoherent mixture of the Zeeman sublevels of the ground 
state. 


It is surprising that, unlike the transmittivity |T^/°^| 2 
for the propagation direction /o, the transmittivity 
|T^ /o) | 2 for the propagation direction — f 0 does not re¬ 
duce to zero with increasing atom number N. Due 
to this feature, the atomic array can act as an optical 
diode even in the limit IV —> oo [34]. The property 

|| — > | T ^* 0 ' 1 1 ^ 0 in the limit N —> oo is not related 
to the EIT effect. Indeed, this property occurs even in the 
case where there is no control field (fl c = 0). This prop¬ 
erty is a consequence of the Bragg resonance condition 
and the difference between the coupling coefficients |G+| 
and |G_| for the guided modes with the opposite prop¬ 
agation directions / = + and / = —, respectively. The 
difference between |G+| and |G_| is related to the exis¬ 
tence of the l ong itudinal component e z of the x-polarized 
guided field [55]. 

In order to get insight into the above result, we con¬ 
sider the change from the case of N atoms to the case 
of IV + 1 atoms by adding an atom to the array of N 
atoms under the Bragg resonance condition /3 p A = nir. 
We use the recurrence relations «5ZD to describe the 
change. In the limit N —> oo, we have Rn —> Roo and 
0. We find that a nonzero asymptotic solution 

■> (-1 )(A+ 1 )n r ^/° ) ^ o exigts if 


T, 


(/o) 


T; 


N 
C-/o) 


N 


jV-io) 

1 — RooR 


(63) 


3. Limiting values 


We derive analytical expressions for the transmittivity 
and reflectivity of the atomic array in the limit of large 
N. 

We first analyze the case of the cc-polarization scheme 
under the Bragg resonance condition. In this case, in the 
limit N —i 00, we obtain the reflection coefficient 


Rn —i Roo — 


M - \e z \ 

\ e r \ + l e z| 


(61) 


and the transmission coefficients 


iii /o) i^o, 
\T { N ~ fo) \ |rj-A>| 


4|e r |M 

(M + M ) 2 


(62a) 

(62b) 


where /o = sign (M e — M g ). 

It is clear that the limiting value R a0 is determined by 
the guided-mode profile functions e r and e 2 only. Since 


that is, if Roo = (1 — ^)/R } in agreement with 

Eq. (1BT1) . The inverse of the denominator 1 — RooR hr 
the expression on the left-hand side of condition (IMl) 
characterizes the enhancement due to multiple reflections 
between the initial array and the added atom. Condi¬ 
tion (1631) requires that the enhancement of transmission 
caused by multiple reflections compensates the reduction 
caused by a single pass through the added atom. Thus, 
in the limit of large IV, the reflection coefficient Rn of the 
array approaches an appropriate value Roo, which satis¬ 
fies condition m , and hence the transmittivity |T^ ^ | 2 
for the propagation direction —f 0 = —sign (M e — M g ) 
tends to a nonzero value given by Eq. (I62bl) . 

We now analyze the case of the y-polarization scheme 
under the Bragg resonance condition. In this case, we 
find 


Rn = 


NR 


1-(IV-1)1? 

T n = (-lJ^+D 


1-(IV-1)1?' 


( 64 ) 
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In the limit TV —> oo, we have 

Rn —> — 1, (65a) 

T n -> 0, (65b) 

that is, |i?7v| 2 —>• 1 and \Tn\ 2 —0. This result means 
that the atomic array under the Bragg resonance con¬ 
dition can act as a perfect mirror for the y-polarized 
guided light fields in the limit TV —> oo. The loss 
due to the scattering into the radiation modes is sup¬ 
pressed due to the collective enhancement of scattering 
into the backward guided modes. In the limit NR -C 1, 
Eqs. (1M1) yield Rn — NR + TV (TV — 1 )R 2 and Tn — 
(_l)(JV+i)«[T + (TV - T)RT + (TV - l) 2 i? 2 T], The last 
terms in these expressions contain TV 2 . They are signa¬ 
tures of the collective effects. 


4- Dependence on the array period 

We consider the effect of the array period A on the 
Bragg resonance. For this purpose, we plot in Figs. EU1 
andl2Tlthe dependences of the transmittivity |X^/ | 2 , the 

reflectivity |T?at| 2 , the transmitted-field group delay t?\ 
and the reflected-field group delay tr on the array period 

A. 



Array period A (nm) 


FIG. 20. Dependences of the transmittivities |T]y^| 2 and 
|T)^ ^| 2 , the reflectivity \Rn\ 2 , and the group delays 

and tr on the array period A in the case of the x- 
polarization scheme. The detuning of the probe held is A = 0. 
The number of atoms in the array is TV = 200. Other param¬ 
eters are as for Fig. [2] 

The figures show that both the transmittivity and the 
reflectivity have a local maximum at the array period 
A = 745.16 nm, which satisfies the Bragg resonance con¬ 
dition for the field frequency lo p = ojq- The coexistence of 
the local maxima of the transmittivity and reflectivity at 
the Bragg resonance is an interesting feature. This result 
indicates that the scattering from the atoms into the ra¬ 
diation modes is suppressed due to the Bragg resonance, 



Array period A (nm) 


FIG. 21. Dependences of the transmittivity |TW| 2 , the re¬ 
flectivity | Rn | 2 , and the group delays tt and tr on the array 
period A in the case of the {/-polarization scheme. The detun¬ 
ing of the probe field is A = 0. The number of atoms in the 
array is TV = 200. Other parameters are as for Fig. [3] 

in agreement with the results of Ref. [51] . It is clear that, 
in the vicinity of the Bragg resonance, the reflectivity for 
the //-polarized guided field [see Fig. l2TT c')] is larger than 
that for the z-polarized guided field [see Fig. l20Te)]. We 
note that the linewidth of the dependence of the reflec¬ 
tivity and transmittivity on the array period A is on the 
order of a few nanometers. The reason is that, when 
f3 p 5A -C 1, that is, when SA -C Xp = 2tt//3 p — 745.16 
nm for X p = Xq = 852.35 nm, there is no significant de¬ 
viation from the Bragg resonance condition /3 p A = nir , 
with n = 1,2,.... Here, SA is a small deviation of the 
array period A from a Bragg resonant value. We observe 
narrow dips in the A dependence of tr [see Figs. EOT f'l 
and r2TT cl 1 ]. These dips occur at the points where |i?jv| 2 
have minima [see Figs. 0e) and [2lT cl]. In other words, 
the dips correspond to the interference fringes that sur¬ 
round the Bragg resonance. The separation AA between 
the positions of the dips is determined by the half-period 
of the function sinTV£, which appears in Eqs. ([561). Us¬ 
ing the approximation £ — /3 P A, we obtain the estimate 
A A ~ n/N/3 p . 

5. Dependence on the control field intensity 

In order to show the effect of the control field on the 
transmission and reflection of the probe field, we plot in 
Figs. [HI and [55] the transmittivity |Tj^| 2 , the reflectivity 

| Rn ] 2 , the transmitted-field group delay Tjf\ and the 
reflected-field group delay tr as functions of the intensity 
I c of the control field £ c . 

The left columns of these figures show that, when I c 
increases, the transmittivity |T)^' ) | 2 increases but the re¬ 
flectivity |T?at | 2 decreases. We observe from the right 
columns of Figs. [?21 and E51 that, when I c increases and 
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Control field intensity I (mW/cm 2 ) 


FIG. 22. Dependences of the transmittivities IT^'l 2 and 
|T^ J | 2 , the reflectivity \Rn\ 2 , and the group delays 

and tr on the control held intensity I c in the case of the 
axpolarization scheme. The period of the array is A = 745.16 
nm, which satisfies the second-order Bragg resonance condi¬ 
tion. The detuning of the probe held is A = 0. The number 
of atoms in the array is N = 200. Other parameters are as 
for Fig. |2] 



Control field intensity / (mW/cm 2 ) 

FIG. 23. Dependences of the transmittivity |Tiv| 2 . the re¬ 
flectivity |i?jv| 2 , and the group delays tt and tr on the con¬ 
trol held intensity I c in the case of the y-polarization scheme. 
The period of the array is A = 745.16 nm, which satishes the 
second-order Bragg resonance condition. The detuning of the 
probe held is A = 0. The number of atoms in the array is 
N = 200. Other parameters are as for Fig. [3] 


( f) 

is not too small, the group delay Tfl of the transmitted 
field decreases, and the absolute value |rp| of the nega¬ 
tive group delay tr < 0 of the reflected field increases. 
The above-mentioned features are the results of the dis¬ 
persion properties in the EIT window. 


6 . Slow transmitted light and fast reflected light 

We plot in Figs. [24] and [25] the time dependences of 
the normalized intensities |-4out,ref/-4o| 2 of the trans¬ 
mitted/reflected pulses in the cases of the x- and y- 
polarization schemes, respectively. The parameters used 
are as for Figs. [T5] and [13] except for the array period 
A = 745.16 nm, which satisfies the second-order Bragg 
resonance condition. 



Time t (ps) 


FIG. 24. (Color online) Time dependences of the normal¬ 
ized intensities |-4 ou t,ref/-4o| 2 of the transmitted and reflected 
pulses in the case of the axpolarization scheme. The input 
pulse is of the Gaussian form, with the central frequency 
uj p = wo, the pulse length to = 2 ps, and the peak time 
t = 0. The number of atoms is N = 200. The array period 
is A = 745.16 nm, which satisfies the second-order Bragg res¬ 
onance condition. Other parameters used are as for Fig. [2] 
The dotted red line is for the input-pulse peak time t = 0 and 
is a guide to the eye. 

We observe that, when the probe pulse is x polarized 
and incident onto the array in the direction f = + [see 
Fig. f23T a)] and when the probe pulse is y polarized an 
incident in an arbitrary propagation direction / = ± [see 
Fig. the transmitted pulse is significantly weak¬ 

ened and delayed by the atomic array. For N = 200 
atoms in the array with the length L = (N — 1)A ~ 148 
pm, we obtain the group delays of about 1.11 ps in the 
case of Fig. VIM a) and 0.37 ps in the case of Fig. l25T a). 
These group delays correspond to the group-velocity re¬ 
duction factors c/Vg ~ 2.2 x 10 6 and ~ 0.75 x 10 6 in 
the cases of Figs. !23T a) and I25f a). respectively. Figure 
123T b) shows that, when the probe pulse is x polarized 
and incident onto the array in the direction / = —, the 
intensity reduction and the group delay of the transmit¬ 
ted pulse are not significant. According to Figs. l23T c) 
and 125T b). the reflected pulse has a negative group de¬ 
lay (fast light) and its intensity is substantial due to the 
Bragg resonance condition. 

Comparison between Figs. [23] and [32] and between 
Figs. 1251 and fl3l shows that the shapes and group delays 
of the transmitted and reflected pulses are not affected 
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Time t (ps) 


FIG. 25. (Color online) Time dependences of the normal¬ 
ized intensities |-4 0 ut,ref/-4o| 2 of the transmitted and reflected 
pulses in the case of the j/-polarization scheme. The input 
pulse is of the Gaussian form, with the central frequency 
uip = wo, the pulse length tq = 2 fis, and the peak time 
t = 0. The number of atoms is IV = 200. The array period 
is A = 745.16 nm, which satisfies the second-order Bragg res¬ 
onance condition. Other parameters used are as for Fig. [2] 
The dotted red line is for the input-pulse peak time t = 0 and 
is a guide to the eye. 

much by the Bragg resonance condition. However, the 
intensity of the reflected pulse in the presence of a Bragg 
resonance [see Figs. l24T c) and l25T b)] is much higher than 
that in the absence of a Bragg resonance [see Figs. Me) 
and lTSI b)]. 

E. Photonic band gaps 

It is known that, in the neighborhood of a Bragg res¬ 
onance, where A = A res = n\p/2 with n = 1,2,..., 
band gaps may be formed when the number of atoms 
N in the array is large enough HH, |45| [51 ]. In order to 
show the band gaps, we plot in Figs. [26] and [27] the trans- 
mittivity and the reflectivity |i?Ar| 2 of the atomic 

array with the period A = 745.16 nm (corresponding to 
A A = A — A res = 0 for A p = Ao) as functions of the field 
detuning A for a very large number of atoms, namely 
for N = 200 000. In the calculations for these figures, 
we used a large value for the control field intensity 7 C , 
namely I c = 20 mW/crn 2 , in order to show clearly the 
tiny features of the transmittivity and reflectivity for the 
probe field in the EIT window. Other parameters are as 
for Figs. [2] and [3] 

Figure Uni shows the frequency dependences of |T^ f ' l | 2 , 
|T^| 2 , and \Rn\ 2 for the x-polarization scheme in the 
vicinity of a Bragg resonance A A = 0. We observe from 
this figure that, in addition to a narrow plateau around 
the atomic resonance A = 0, there are two wide plateaus, 
one on the left side and the other one on the right side. 
The left- and right-side plateaus are the photonic band 




FIG. 26. Photonic band gaps in the case of the ^-polarization 
scheme with AA = 0. The transmittivities |T^, +) | 2 (a) and 
\T ( n ^| 2 (b) and the reflectivity |i?jv | 2 (c) are plotted as func¬ 
tions of the detuning A. The period of the array is A = 745.16 
nm, which satisfies the second-order Bragg resonance condi¬ 
tion. The number of atoms in the array is N = 200 000. 
The intensity of the control field is I c = 20 mW/cm 2 (the 
corresponding Rabi frequency is = 2 .O 670 ). Other param¬ 
eters are as for Fig. [2] The insets show the narrow structures 
around the point A = 0. 



Probe field detuning A/2jt (MHz) 


FIG. 27. Photonic band gaps in the case of the y-polarization 
scheme with AA = 0. The transmittivity |Tjv | 2 (a) and the 
reflectivity |Aiv| 2 (b) are plotted as functions of the detuning 
A. The period of the array is A = 745.16 nm, which satisfies 
the second-order Bragg resonance condition. The number of 
atoms in the array is N = 200 000. The intensity of the control 
field is I c = 20 mW/cm 2 (the corresponding Rabi frequency 
is = I. 9470 ). Other parameters are as for Fig. [3] The 
insets show the narrow structures around the point A = 0. 


gaps that extend over the frequency range from u c — A max 
to u c A min and from uj c T A m j n to u c T A max [39l , [5l| . 
Here, we have introduced the notations u c = (u o+wi a t)/2 
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and 


^max — 

A m jn — 


Alt , 

u 0 \e r \ 2 v g 

4 + 

A : 

'*L , 

u 0 \e z \ 2 v g 


A 


( 66 ) 


with <Si a t = wiat - w 0 , and u 0 = uj p dl g /e 0 hv g . The Bragg 
resonant frequency wi a t is determined by the equation 
/3(wi a t)A = mr, where the integer number n = 1,2 ,... is 
the order of the Bragg resonance. In deriving Eqs. (EH) 
for the edges of the band gaps, we have neglected the 
atomic decay rate T eg and the control field Rabi fre¬ 
quency Q c . In the case of |<5i a t| <C \/u{)\e r ^ z \ 2 v g / A, we 
find that the band gaps will be formed when N 7V gap , 
where iV gap = ^v g /uoA/(\e r \ — |e 2 |). We obtain from 
Eqs. (liEl) the estimates A min = 27r x 847 MHz and 
A max = 27r x 1544 MHz for the parameters of Fig. [211 
In the band gap regions of Fig. [26l we have |-Rjv | 2 — 1 
and IT^I 2 ~ 0. In the central plateau, we have 
|R W | 2 ~ 0.085 < 1, |T ^ +) | 2 ~ 0, and |T ^ _) | 2 ~ 0.84, in 
agreement with Eqs. fl§H) and (l62l) . The limiting values 
jTic^l 2 ~ 0 and |Tio ^| 2 ~ 0.84 ^ 0 of the transmittivi- 
ties in the central plateau area indicate that the atomic 
array can operate as an optical diode even in the limit of 
an infinitely large value of N. The insets of Figs. 1251' a’)- 
EH' c'l show that there are narrow structures in the fre¬ 
quency dependences of the transnrittivity \T^\ 2 and the 
reflectivity |Rat | 2 in the vicinity of the atomic resonance 
A = 0. These tiny structures are the consequences of the 
interplay between the scattering into the guided and ra¬ 
diation modes, the EIT effect, and the Bragg resonance. 

Figure [271 shows the frequency dependences of |T/v | 2 
and |-Rat | 2 for the y-polarization scheme in the vicinity 
of a Bragg resonance A A = 0. We observe that there is a 
wide plateau around the atomic resonance A = 0. This 
plateau corresponds to the set of two photonic band gaps 
that extend over the frequency range from oj c — A max to 
w c -|-A max , with frequencies between the atomic frequency 
cuq and the frequency wi a t excluded [39, [5lj| . Here, we 
have introduced the notation 


A m „ = + (67) 

In deriving Eqs. E3 for the edges of the band gaps, 
we have neglected the atomic decay rate T eg and the 
control field Rabi frequency fi c . In the case of |<5i at | -C 
y/ u o\ e ip\ 2 v g / A, we find that the band gaps will be formed 
when N » IV gap , where A gap = ^2v g /3u 0 A/\e (p \. We 
obtain from Eq. (1671) the estimate A max = 27r x 1561 
MHz for the parameters of Fig. [27] In the band gap re¬ 
gion of Figs. [271 we have |Rjv | 2 ~ 1 and |T/v | 2 — 0, in 
agreement with Eqs. (1651) . The narrow structures in the 
insets of Figs. [771 a) and lTTI bl are the consequences of the 
interplay between the scattering into the guided and ra¬ 
diation modes, the EIT effect, and the Bragg resonance. 
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FIG. 28. As Fig. 1261 but for AA = 0.3 nm. 



Probe field detuning A/2jt (MHz) 


FIG. 29. As Fig. 1271 but for AA = 0.3 nm. 


The band gaps illustrated in Figs. [7711 and [771 may ap¬ 
pear even when Q c = 0 , that is, when there is no control 
field j3§, 0, [Elj. These band gaps are not related to the 
EIT effect. We call them the non-EIT band gaps. When 
the condition 


I r | ^ / tfco \&r ip zl^Vg 

IM > Y ^ »7o (68) 

is satisfied, one of the two non-EIT band gaps is far away 
from the atomic resonance frequency uiq while the other 
one is close to ujq. We plot in Figs. [751 and [751 the trans- 
mittivity |T) 1/' > | 2 and the reflectivity |R(v | 2 of the atomic 
array as functions of the field detuning A in the case of 
A = 745.46 nm, that is, AA = A — A res = 0.3 nm, which 
corresponds to <5i at = —v g /3 2 AA/nn ~ —107 GHz. Other 
parameters are the same as for Figs. l26l and 1271 We ob¬ 
serve from the insets of Figs. E5T c') and EHb) that, in 
the vicinity of A = 0, there are two peak regions where 
the reflectivity |Rat | 2 is significant [43, (48(. One of these 
peak regions is broad and is a non-EIT band gap. The 
other peak region is narrow and is also a band gap. The 
occurrence of this additional band gap is due to the EIT 
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effect caused by the action of the control field £ c [47l.l48j. 

In the above calculations, we did not include the ex¬ 
perimental limited filling ratio of the atomic array. In 
the Vienna atom trap experiment (35j . due to the small 
trapping volumes,the loading operated in the collisional 
blockade regime |61j. This results in an occupancy of at 
most one atom per trapping site. For the parameters of 
the experiment [351, the filling ratio is about 0.5 |35i. [6lj. 
We note that, in the case where the array period A is 
far from the Bragg resonance, the inclusion of the filling 
ratio will lead to just a reduction of the atom-number 
density ua and an increase of the medium length L for 
a given atom number N. In the case where the array 
period A is near to a Bragg resonance, the presence of 
a void in the array will lead to just an additional phase 
of nir for the transmission coefficient Tjj ', where n is 
the order of the Bragg resonance. Therefore, we expect 
that the inclusion of the filling ratio will not affect the 
results for the above limiting cases. In the intermediate 
case where the array period is not far from and not close 
to a Bragg resonance, the inclusion of the filling ratio 
will lead to a reduction of the reflectivity (and possibly 
also a reduction of the transmittivity) of the array with 
a given atom number N. In particular, the widths of the 
resonances in Figs. EU1 and Ell will be reduced. 


V. SUMMARY 

We have studied the propagation of guided light along 
an array of three-level atoms trapped in the vicinity of 
an optical nanofiber under the EIT condition. We have 
examined two schemes of atomic levels and field polariza¬ 
tions where the guided probe field is quasilinearly polar¬ 
ized along the major principal axis x or the minor princi¬ 
pal axis y. We have derived the coupled-mode propaga¬ 
tion equations, the input-output relation, the scattering 
matrix, and the transfer matrix for the transmitted and 
reflected fields. We have taken into account the com¬ 
plexity of the polarization of the guided field and the 


discreteness and periodicity of the atomic positions in 
the array. We have calculated the reflection and trans¬ 
mission coefficients. We have found that, when the array 
period is far from the Bragg resonance, the reflection is 
negligible and the homogeneous-medium approximation 
can be used. We have numerically demonstrated that 200 
cesium atoms in a linear array with a length of 100 //rii 
at a distance of 200 nrn from the surface of a nanofiber 
with a radius of 250 nm can slow down the speed of 
guided probe light by a factor of about 3.5 x 10 1 2 3 4 * 6 (the 
group delay is about 1.17 ys). In the neighborhood of 
the Bragg resonance, a significant fraction of the guided 
probe light can be reflected back with a negative group 
delay (that is, with a positive group advance). The re¬ 
flectivity and the group delay of the reflected field do 
not depend on the propagation direction of the probe 
field. However, when the input guided light is quasi¬ 
linearly polarized along the major principal axis x, the 
transmittivity and the group delay of the transmitted 
field substantially depend on the propagation direction 
of the probe field. When the input guided light is quasi¬ 
linearly polarized along the major principal axis x and 
propagates in the direction / = — sign(M e — M g ), un¬ 
der the Bragg resonance condition for an array of atoms 
in an appropriate internal state, the transmission of the 
guided light is not zero even in the limit of large atom 
number N. This result indicates that the atomic array 
can operate as an optical diode^ even in the limit of in¬ 
finitely large atom numbers [34l ]. The directionality of 
transmission of guided light through the atomic array is 
a consequence of the existence of a longitudinal compo¬ 
nent of the guided light field as well as the ellipticity of 
both the field polarization and the atomic dipole vector. 
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